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Many new cc or charmonium-likc states have been 
discovered at the _B-factories in the energy region 
above DD threshold. Of these, the X(3872) Q], 
X c2 (2F)(3930) [2], V(3940) [3], and V(4260) [4] reso- 
nances are now well-established. Since the V(4260) is 
produced via initial-state radiation (ISR) in the reaction 
e + e~ — > 7isr J/ipir + ir~, it has J PC — 1 . In addition 
to the V(4260), two more J PC = 1 — states, the V(4360) 
and the V(4660), have been reported in ISR production, 
via e+e- -> 7lSR </'(2S)7r + 7r- |6]. The F(4660) has 
been observed only in the Belle experiment [5] , and so it 
is important to confirm the existence of this state. 

In this paper we utilize the ISR mechanism to study 
the reaction e + e~ — > vb(2S)n + n- in the center-of-mass 
(cm.) energy (E cm ) range 3.95 - 5.95 GeV, where the 



ip(2S) decays to J/ipir + Tr~ or to l + l~, with repre- 
senting either e + e~ or 

We use a data sample corresponding to an integrated 
luminosity of 520 fb^ 1 , recorded by the BABAR detector at 
the SLAC PEP-II asymmetric-energy e + e~ collider op- 
erating at and near the cm. energies of the T(nS) (n = 2, 
3, 4) resonances. The detector is described in detail else- 
where [7]. Charged-particle momenta are measured in a 
tracking system consisting of a five-layer, double-sided, 
silicon vertex-tracker (SVT) and a 40-layer central drift 
chamber (DCH), both coaxial with the 1.5-T magnetic 
field of a superconducting solenoid. An internally reflect- 
ing ring-imaging Cherenkov detector, and specific ion- 
ization measurements from the SVT and DCH, provide 
charged-particle identification (PID). A CsI(Tl) electro- 
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magnetic calorimeter (EMC) detects and identifies pho- 
tons and electrons. Muons are identified using informa- 
tion from the instrumented flux-return system. 

We reconstruct events corresponding to the reaction 
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FIG. 1. (a) The ^>(2S)tt + -k~ invariant mass distribution 
from the kinematic threshold to 5.95 GeV/c 2 for ip(2S) — > 
J I ij)n + iv~ ; the points with error bars represent the data in 
the ip(2S) signal region, and the shaded histogram is the 
background estimated from the ip(2S) sideband regions. The 
solid curve shows the result of the fit described in the text. 
The dashed (dotted) curves indicate the individual resonant 
contributions for constructive (destructive) interference, (b) 
The corresponding distributions for ip{2S) — > l + l~. (c) 
The combined ip(2S)ir + -K~ invariant mass distribution for 
ip(2S) -> J/<ipTT + Tv~ and ip(2S) l + l~. The solid curve 
shows the result of the fit. The dashed curve represents the 
background, while the dotted curves indicate the individual 
resonant contributions. There is only one solution in this case. 

e + e~ — > 7isR?/'(2>S')7r + 7r~, where 7isr represents a pho- 



ton that is radiated from the initial-state e^- , thus lower- 
ing the cm. energy of the e + e~ collision which produces 
the ijj(2S)iT + 7r~ system. We do not require observation 
of the ISR photon, since it would be detectable in the 
EMC for only - 15% of the events. 

For the 4>(2S) — > J/ipir + ir~ decay mode, we select 
events containing exactly six charged-particle tracks, and 
reconstruct J/tp candidates via their decay to e + e~ or 
. For each mode, at least one of the leptons must be 
identified on the basis of PID information. When possi- 
ble, electron candidates are combined with photons to re- 
cover bremsstrahlung energy loss in order to improve the 
J/tjj momentum measurement. An e + e~ pair with invari- 
ant mass within (—60,-1-45) MeV/c 2 of the nominal J/ip 
mass [5] is accepted as a J/ip candidate, as is a pair 
with mass within (—45, +45) MeV/c 2 of this value. Each 
J/0 candidate is subjected to a geometric fit in which 
the decay vertex is constrained to the e + e _ collision axis 
within the interaction region; the % 2 -probability of the fit 
must be greater than 0.001. An accepted J/tp candidate 
is kinematically constrained to the nominal J/ip mass [5] , 
and combined with a pion pair to form a J/ipTr + n~ can- 
didate. The J j t/jit + tt~ combinations with invariant mass 
within 10 MeV/c 2 of the nominal ip{2S) mass are taken 
as ip{2S) candidates, and hereafter we refer to this as 
"the ip(2S) signal region". The ^p{2S) candidate is re- 
fit requiring that the x 2 -probability for the vertex fit be 
greater than 0.001. It is then combined with two addi- 
tional pions of opposite charge, each of which is identified 
using PID information, to reconstruct a ?p(2S)Tr + Tt~ can- 
didate. A further geometric fit with the ip(2S) candidate 
mass-constrained to the nominal mass value j5] is per- 
formed. Candidates with x 2 -fit probability greater than 
0.001 are retained for further analysis. 

For the decay mode ip(2S) l + l~, we select events 
containing exactly four charged-particle tracks, and re- 
construct ip(2S) candidates via their decay to e + e~ or 
fi + fi~ . An e + e~ (/i + /i~) pair with invariant mass within 
(-40,+30) MeV/c 2 ((-30,+30) MeV/c 2 ) of the nominal 
ijj(2S) mass [8] is accepted as being within the ip(2S) sig- 
nal region. Each such candidate is subjected to the same 
geometrical fit and mass constraint procedure as applied 
for the ip{2S) — > J/^7r + 7r~ mode. A surviving candidate 
is combined with a pion pair to form a tJj(2S)tt + tt~ can- 
didate. 

For ip{2S) -> J/^tt+tt- (4>(2S) -> l+r), the dif- 
ference between the cm. momentum of the hadronic 
iP(2S)tt + it~ system and the value expected for an ISR 
event (i.e. (s — m 2 )/2y / s, where m is the ip(2S)ir + ir~ in- 
variant mass) must be in the range (—0.10,-1-0.70) GeV/c 
((-0.70,+0.60) GeV/c) to be consistent with an ISR pho- 
ton. We require the transverse component of the missing 
momentum to be less than 2.0GeV/c (1.7GeV/c). If the 
ISR photon is detected in the EMC, its momentum vector 
is added to that of the i/j(2S)tt + it~ system in calculat- 
ing the missing momentum. For the events for which 
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TABLE I. Results of the fit to the ifj(2S)Tr + -K~ invariant mass 
distributions for ip(2S) — ¥ J/i}TT + ir~ . The first errors are 
statistical and the second systematic; B x r ee is the product 
of the branching fraction to ^(2S')7r + 7r _ and the e + e~ partial 
width (in eV), and <f> is the relative phase between the two 
resonances (in degrees). 

Parameters First Solution Second Solution 

(constructive (destructive 
interference) interference) 
4340 ± 16 ± 9 
94 ± 32 ± 13 
6.0 ± 1.0 ± 0.5 7.2 ± 1.0 ± 0.6 
4669 ± 21 ± 3 
104 ± 48 ± 10 
2.7 ±1.3 ±0.5 7.5 ±1.7 ±0.7 
12 ± 27 ± 4 -78 ± 12 ± 3 



Mass Y(4360) (MeV/c 2 ) 
Width Y(4360) (MeV) 
B x r ee (y(4360)) (eV) 

Mass Y(4660) (MeV/c 2 ) 
Width Y(4660) (MeV) 
B x r ee (y(4660)) (eV) 
<M°) 



TABLE II. Results of the fit to the combined ^(2S)tt + tt~ 
invariant mass distributions for ip(2S) — ¥ J/ij)ix + -K^ and 
"4>{2S) — ¥ l + l~ . The first errors are statistical and the second 
systematic; B x F ee is the product of the branching fraction 
to tp(2S)'K + 'K~ and the e + e~ partial width (in eV), and (f) is 
the relative phase between the two resonances (in degrees). 



Parameters 
Mass Y(4360) (MeV/c 2 ) 
Width Y(4360) (MeV) 
B x r ee (Y(4360)) (eV) 
Mass Y(4660)(MeV/c 2 ) 
Width Y(4660) (MeV) 
B x r ee (Y(4660)) (eV) 
4>{°) 



Solution 



4318^i'a± 



-19 

20 ± 



123 ± 
7.4 ± 0.9 

7+6 



3 

13 

0.7 



36l? 2 ±4 



1.4 ± 

25 ± 



14 

0.5 
21 



± 0.2 

± 2 



ip(2S) — ¥ e + e~ , the candidate ir + ir~ system has a small 
contamination due to e + e~ pairs from photon conver- 
sions. We compute the pair invariant mass m e + e - , with 
the electron mass assigned to each pion candidate, and 
remove candidates with m e + e - < 100 MeV/ c 2 . 

For events with multiple ^p{2S) candidates, we select 
the combination that has candidate mass closest to the 
ijj(2S) nominal mass value [5]. We estimate the remain- 
ing background for the decay mode ijj(2S) — ¥ J/ipn + ir~ 
using events that have a J/-07t + 7t~ mass in either of the 
V>(2S) sideband regions (3.566, 3.666) or (3.706, 3.806) 
GeV/c 2 . For the decay mode ip(2S) —¥ e + e~, the cor- 
responding regions are (3.476, 3.576) and (3.776, 3.876) 
GeV/c 2 , while for ip(2S) —¥ the sideband regions 

are (3.516, 3.596) and (3.776, 3.856) GeV/c 2 . 

Figure [l] shows the ^(2S)ir + u~ invariant mass distri- 
butions for the selected ip(2S) events corresponding to 
the decays (a) ip(2S) ->■ J/^tt+tt-, (b) ip(2S) ~¥ l+l~, 
and (c) the combined sample for ip(2S) — ¥ J/ipn + TT~ 
and ip(2S) — ¥ . The background is estimated from 
the ip(2S) mass sidebands as described above. In Fig. 
[ljtwo structures are evident, the first near 4.35 GeV/c 2 , 
and the second near 4.65 GeV/c 2 . We attribute these 
peaks to the Y(4360) [5] and to the F(4660) [6J, respec- 



tively. We first fit the distribution shown in Fig. [ija) in 
order to extract the parameter values of the resonances. 
We then perform a second fit to the combined J/07r + 7r~ 
and data of Fig. |ljc), where the signal yields are 

larger, but where these come at the cost of the large 
background associated with the dilepton channels. For 
both distributions we perform an unbinned, extended- 
maximum-likclihood fit to the ^>{2S)it + it~ mass distri- 
bution from the signal region, and simultaneously to the 
background mass distribution. We describe the latter 
by a fourth-order polynomial in tJ){2S)tt + mass, m, 
for the fit to the data of Fig. [IJa), and by a third-order 
polynomial for the fit to the data shown in Fig. [ijc) . 

The mass dependence of the signal function is given by 
/(m) = e(m) • C(m) ■ cr(m); e(m) is the mass-dependent 
signal-selection efficiency obtained from a MC simula- 
tion which uses a ip(2S)ir + Tr~ phase space distribution; 
its value increases from 1% at 3.95 GeV/c 2 to 12% at 
5.95 GeV/c 2 for ip(2S) ->■ J/^tt+tt" and from 1% at 
3.95 GeV/c 2 to 14% at 5.95 GeV/c 2 for ip(2S) -> 
The function £(m) is the mass-distributed luminosity [Hj 
(we ignore the small corrections due to initial-state emis- 
sion of additional soft photons); £(m) increases from 102 
pb-750 MeV to 202 pb"V50 MeV from 3.95 GeV/c 2 to 
5.95 GeV/c 2 . 

The cross section, <r(m), is described by the follow- 
ing function, which takes into account the possibility of 
interference between the two resonant amplitudes, since 
they have the same quantum numbers (J 



pc — i — 



I"): 



a ( m ) = 1 ^L.\A 1 (BW) 



A 2 (BW) 




(1) 



PS(m 2 ) 



where C = 0.3894 • 10 9 GeV 2 pb, and PS(m) represents 
the mass dependence of ip(2S)ir + ir~ phase space; 4> is the 
relative phase between the amplitudes A\ and A 2 . The 
complex amplitude Aj is written as 



A 3 {BW) = 



r 



i/)(2S)ir+7r-)j 



irrijTj 



(2) 



where nij is the resonance mass and T, its total width; 



(r e + e - -r 



is the product of the partial widths 



to e + e~ and to -0(2S')7r + 7r _ . 

In the fit procedure f(m) is convolved with a Gaus- 
sian resolution function obtained from MC simulation. 
This function has root-mean-squared (r.m.s.) deviation 
which increases linearly from 2MeV/c 2 at ~ 3.95 GeV/c 2 
to 5 MeV/ c 2 at - 5.95 GeV/c 2 . In the likelihood function, 
when the fit is performed to the ip(2S) — ¥ J/?pTr + ir~ data, 
ct(to) is multiplied by B{ip{2S) -» J /^tt+tt-) x B(J/ip 

since the fit is to the corresponding observed 
events. Similarly, for ip(2S) — ¥ l + l~, a(m) is multi- 
plied by B(ij)(2S) —¥ where I = e or /1, in fitting 
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FIG. 2. (a) The comparison between the observed 
■4}(2S)it + it~ (ip(2S) — > J/tf}TT + n~) invariant mass spectrum 
from BABAR (dots) and that from Belle (hatched histogram), 
(b) The combined BABAR and Belle ip{2S)-K + it~ (ip(2S) 
J/i/m + iT~) invariant mass spectrum. 
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FIG. 3. The cross section for the reaction e + e~ — > 
ip(2S)ir + -K~ as a function of cm. energy obtained by us- 
ing Eq. (|3| (points with error bars); the curve shows the cm. 
energy dependence which results from the fit to the data of 
Fig. ga). 



the data. The results of the fits are shown in Fig. [TJa) 
and in Fig. [ljc), and the extracted parameters are sum- 
marized in Tables [i] and [TTJ respectively. The signifi- 
cance of the F(4660) signal for both fits is greater than 
5cr where a is the standard deviation. For the fit to 
the distribution in Fig. [lja), we obtain two solutions, 



one corresponding to constructive interference and one 
to destructive interference between the resonant ampli- 
tudes. The mass and the width values of the resonances 
are the same for each solution. However, the values of 
r e + e - x B(ip(2S) — > J/y>7T + 7r~) and (j) are different (see 
Table [i]), although the maximum likelihood value is ex- 
actly the same for each fit. For the fit to the distribution 
in Fig. [ljc) , only a solution showing constructive inter- 
ference is obtained, for which the parameter values are 
consistent within error with those for the first solution 
in Table [I] The results summarized in Table [I] agree well 
with those obtained in the Belle analysis [B], for which 
the data sample is about the same size as that for the 
i/j(2S) — > J/V>7T + 7r _ decay mode in the present analysis 
(see Fig.[2ja)). We infer that, even if our data sample for 
this mode were doubled in size, the ambiguity in the fit 
results would persist. The inclusion of the ip(2S) dilep- 
ton decay modes increases the signal sample by 40% over 
that for the ip(2S) — > J/ipir + ir~ mode alone. This in- 
crease is obtained at the cost of introducing a background 
contribution which is larger by 50% than the combined 
signal sample. The fit to this sample yields only one so- 
lution (Table [Tl| . However, the comparison of our results 
for the ip{2S) — > J/ipir + ir~ analysis to those from the 
Belle analysis [6J leads us to conclude that the apparent 
resolution of the fit ambiguity is due, not to the slightly 
increased signal sample, but rather to the presence of the 
large background shown in Fig. [ljc). For this reason wc 
discount the results summarized in Table [Hj and confine 
our attention to the results from ip(2S) — > J/i/jir + ir~ de- 
cay in the remainder of the analysis. 

The fit results of Tableland the tP{2S)tt + 7t- invariant 
mass spectrum of Fig[2ja) agree very well with those ob- 
tained by the Belle Collaboration [B]. Each distribution 
(Fig. [5Ja)) shows evidence of two resonant signals (note 
that the Belle distribution ends at 5.5 GeV/c 2 ). This is 
even more apparent in Fig.[2jb), where we have added the 
distributions to obtain a mass spectrum corresponding to 
an integrated luminosity of ~ 1.2 ab _1 . The existence of 
two structures is quite clear, and there is even a hint of 
some activity in the vicinity of 5 GeV/c 2 . 

For the decay mode ~> J/4' 1T+1T ~ ^ we calculate 

the e + e~ — > iP(2S)tt + it~ cross section after background 
subtraction for each ■!/;(2S')7T + 7r~ mass interval, i, using 



^obs 



bkg 



e,; ■ Li-B 



(3) 



where n° hs is the number of observed events, n^ kg is the 
number of background events, q is the average efficiency, 
and Li the integrated luminosity [3] for interval i; B rep- 
resents the product B(ip(2S) -> J/ipir + Tr~) ■ B(J/ip -> 
The resulting dependence of the cross section on 
cm. energy is shown in Fig. [3] We sum over the data 
points in Fig. [3] and obtain a model-independent inte- 
grated cross section value of 311^gg (stat) ±11 (syst) pb 



for the region 3.95-5.95 GeV. The curve shown in Fig. [3] 
results from the fit to the data of Fig. [ija), and provides 
an adequate description of the measured cross section. 



TABLE III. Systematic uncertainty estimates for the param- 
eters used in the fit to the data of Fig. [TJa). 
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FIG. 4. The n + n~ invariant mass spectrum for the ip{2S) — > 
J/ipn + ir~ channel in the ij){2S)-K + -r~ mass region (a) 4.0-4.5 
GeV/c 2 , and (b) 4.5-4.9 GeV/c 2 . The histogram represents 
a MC distribution corresponding to the decay according to 
phase space of (a) one resonance with a mass of 4.360 GeV/c 2 
and width 70 MeV, and (b) one resonance with a mass of 4.660 
GeV/c 2 and width 50 MeV. Each histogram is normalized to 
the corresponding data sample. 

Our estimates of systematic uncertainty result from the 
sources listed in Tables |III| and |IV| where we include the 
latter table, which corresponds to the combined ip(2S) 
decay modes, for the sake of completeness. 

The systematic uncertainties on the fitted values of the 
T(4360) and the Y"(4660) parameters include contribu- 
tions from the fitting procedure (evaluated by changing 
the fit range and the background parametrization), the 
uncertainty in the mass scale (which results from the un- 
certainties associated with the magnetic field and with 
our energy- loss correction procedures [IQlHl]), the mass- 
resolution function, and the change in efficiency when the 
dipion distribution is simulated using the histograms in 
Fig. [4j Uncertainties associated with luminosity, track- 
ing, efficiency and PID affect only T e + e ~ ■ B, and their 
net contribution is 3.3%. Uncertainties on the relevant 
branching fraction values [5] are indicated in Tables 111 
and |IV[ and are relevant only for r e + e - ■ B. These esti- 
mates of systematic uncertainty are combined in quadra- 



Source 


r e+e - -B 


r e+e - • b 


Mass 


r 




(%) 


(%) 


( MeV/c 2 ) 


(MeV) 




(constructive 


(destructive 








interference) 


interference) 






Fit procedure 










for the Y(4360) 


±2.5 


± 1.4 


±9 


± 13 


Fit procedure 










for the Y(4660) 


±14 


± 3.3 


±3 


± 10 


Mass scale 


- 


- 


±0.5 




IVT?!'^ rpenlnfinn 








±1.3 


MC dipion model 


±6.8 


±6.8 






"(*(2S)-»J/V'7r + ir-) 


±1.2 


±1.2 






B(J/4>^l+l~) 


±0.7 


±0.7 






PID, luminosity 










and tracking 


±3.3 


±3.3 






Total (Y(4360)) 


±8 


±8 


±9 


±13 


Total (Y(4660)) 


±16 


±6 


±3 


±10 



TABLE IV. Systematic uncertainty estimates for the param- 
eters used in the fit to data of Fig. |TJc) . 



Source 


r e+e - • B 


Mass 


r 




(%) 


(MeV/c 2 ) 


(MeV) 


Fit procedure 








for the Y(4360) 


±4 


±3 


± 2 


Fit procedure 








for the Y(4660) 


±14 


±2 


± 3 


Mass scale 




±0.5 




Mass resolution 






±1.3 


MC dipion model 


±6.8 






"(V>(2S)->J/»/<7r+T-) 


±1.2 






->>!+;-) 


±0.7 








±1.6 






PID, luminosity 








and tracking 


±3.3 






Total (Y(4360)) 


±9 


±3 


±3 


Total (Y(4660)) 


±16 


±2 


±4 



ture to obtain the values which we quote for the T(4360) 
and F(4660) states. 

In Fig. [4] we show the 7r + 7r~ invariant mass distribu- 
tions for events in the ip(2S)ir + ir~ , i/j(2S) —> J/ipir + ir~ 
invariant mass regions (a) 4.0 GeV/c 2 < m^^gw+n- < 
4.5 GeV/c 2 , and (b) 4.5 GeV/c 2 < m^ 2S ^+ w - < 
4.9 GeV/c 2 . The distributions are consistent with previ- 
ous measurements [6] . In each case, the mass distribution 
appears to differ slightly from the phase-space expecta- 
tion, as shown by the corresponding histogram. For the 
higher mass resonance, there is some indication of an ac- 
cumulation of events in the vicinity of the /o(980) state. 
Similar behavior is observed in [6], and both distribu- 
tions bear a qualitative resemblance to the dipion invari- 
ant mass spectrum from the decay Y"(4260) — ► J/ipir + ir~ 
|12) . The small number of events involved precludes the 
drawing of any definite conclusion. 
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In summary, we have used ISR events to study the 
reaction e + e~ — >• i[)(2S)tt + it~ in the cm. energy range 
3.95-5.95 GeV. We observe two resonant structures, 
which we interpret as the K(4360) and the Y~(4660), re- 
spectively. For the K(4360) we obtain m = 4340 ± 
16 ± 9 MeV/c 2 and T = 94 ±32 ± 13 MeV, and 
for the K(4660) m = 4669 ± 21 ± 3 MeV/c 2 and 
r = 104 ±48 ± 10 MeV; in each case the first uncertainty 
is statistical and the second is systematic. We thus con- 
firm the report in Ref. [H] of a structure near 4.65 GeV/c 2 , 
and obtain consistent parameter values for this state. If 
we include the K(4260), which decays to J/ipir + ir~ [I], 
three charmonium-like states with J PC = 1 have been 
observed in the mass region 4.2-4.7 GeV/c 2 , none of 
which has a well-understood interpretation. 
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